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Abstract 

We address the issue of adapting optical images-based edge detection techniques 
for use in Polarimetric Synthetic Aperture Radar (PolSAR) imagery. We mod¬ 
ify the gravitational edge detection technique (inspired by the Law of Universal 
Gravity) proposed by Lopez-Molina et al, using the non-standard neighbour¬ 
hood configuration proposed by Fu et al, to reduce the speckle noise in polari¬ 
metric SAR imagery. We compare the modified and unmodified versions of the 
gravitational edge detection technique with the well-established one proposed 
by Canny, as well as with a recent multiscale fuzzy-based technique proposed by 
Lopez-Molina et al. We also address the issues of aggregation of gray level im¬ 
ages before and after edge detection and of filtering. All techniques addressed 
here are applied to a mosaic built using class distributions obtained from a 
real scene, as well as to the true PolSAR image; the mosaic results are as¬ 
sessed using Baddeley’s Delta Metric. Our experiments show that modifying 
the gravitational edge detection technique with a non-standard neighbourhood 
configuration produces better results than the original technique, as well as the 
other techniques used for comparison. The experiments show that adapting edge 
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detection methods from Computational Intelligence for use in PolSAR imagery 
is a new field worthy of exploration. 

Keywords: Edge detection, SAR images, Computational Intelligence, 
Gravitational method 


1. Introduction 

Edge detection seeks to identify sharp differences automatically in the infor¬ 
mation associated with adjacent pixels in an image [L. Edge detection for opti¬ 
cal images is nowadays quite an established field. It is traditionally carried out 
using gradient-based techniques, such as the well-known Canny algorithm [2]. 
Techniques based on Computational Intelligence have also been proposed in 
the recent literature. Sun et al [3 proposed the gravitational edge detection 
method, inspired by Newton’s Universal Law of Gravity. Lopez-Molina et al [4] 
proposed a fuzzy extension for this technique, allowing the use of T-norms, a 
large class of fuzzy operators; they also proposed small modifications in the basic 
formalism (see Section [3|. Dankova et al [5] proposed the use of a fuzzy-based 
function, the F-transform; the original universe of functions is transformed into 
a universe of their skeleton models (vectors of F-transform components), mak¬ 
ing further computations easier to perform. Barrenechea et al [6] proposed the 
use of interval-valued fuzzy relations for edge detection, using a T-norm and 
a T-conorm to produce a fuzzy edge image, that is then binarized. This ap¬ 
proach was extended by Chang and Chang |7{. First of all, two new images 
are created—one rather dark and the other rather bright—by applying two dif¬ 
ferent parameters on the linear combinations of the images obtained using min 
and max operators, respectively. Then, the fuzzy edge image is created by 
the difference between these two new images. Another recent approach from 
Computational Intelligence is the multiscale edge detection method proposed 
by Lopez-Molina et al [8|, using Sobel operators for edge extraction and the 
concept of Gaussian scale-space. 

SAR sensors are not as adversely affected by atmospheric conditions and the 
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presence of clouds as optical sensors. Moreover, unlike the optical counterparts, 
SAR sensors can be used at any time of day or night. For these reasons, remote 
sensing applications using SAR imagery have been growing over the years ||. 
SAR images, however, contain a great amount of noise, known as speckle , that 
degrades the visual quality of the images. Caused by inherent characteristics of 
radar technology, this multiplicative non-Gaussian noise is proportional to the 
intensity of the received signal. 

Contrary to what happens with optical images, there are still few algorithms 
specifically dedicated to SAR images m- One interesting means to create edge 
detection algorithms for SAR images is to modify those created for optical im¬ 
ages. However, the use of these methods on SAR images is not straightforward, 
due to speckle. One can either adapt optical image techniques to meet SAR 
data properties, or first preprocess the images using filters and then apply the 
original optical techniques. 

The main purpose of our study is to investigate the application of the grav¬ 
itational edge detection, Here we modify the original 3x3 window: the value in 
each cell in the window is no longer the original one, but the aggregation of a 
set of neighbouring pixels, according to the larger 9x9 neighbourhood configu¬ 
ration proposed by Fu et al m- We propose a typology of experiments to study 
the behaviour of the modified edge detection method, considering polarization, 
image aggregation, and image binarization. We focus on the use of the following 
processes: DAB (edge Detection on non-binary images, Aggregation of the re¬ 
sulting non-binary images, Binarization) and ADB (Aggregation of non-binary 
images, edge Detection on the resulting non-binary image, Binarization). 

We also investigate the use of noise-reduction filters in preprocessing the 
images, by making use of the the well-known Enhanced Lee filter m and a 
filter recently proposed by Torres et al [12]. 

Barreto et al [13] describe a classification experiment, based on a full po- 
larimetric image from an agricultural area in the Amazon region in Brazil. In 
that study, the authors estimated the parameters for probability distributions 
associated to each of the classes of interest, such as water and different types of 
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vegetation and their phenology. They assessed their results in an image formed 
by a mosaic of the classes, with pixel values generated using the parameters 
found for each class. We apply all techniques addressed in this study on twenty 
simulated mosaics, using the parameters estimated in m, considering ampli¬ 
tude images derived from different polarizations. We assess the quality of the 
results, according to Baddeley’s Delta Metric (BDM) [14] . 

We also apply the methods on the real images, but assessment is only visual. 
We compare our results with those produced by the use of Canny’s algorithm [2j 
and the recently proposed multiscale method by Lopez-Molina et al [8]. 

The present study is an extended version of [23 ? in which some of the main 
ideas of this paper were first delineated. However, the present study and m 
differ in the scope of the proposed approach as well as in the reliability of the 
results. Indeed, in HE only one simulated image was used in the experiments 
and only Canny’s technique was compared to its results. Moreover, in the 
previous paper we only addressed the edge detection of the image resulting 
from the aggregation of the three simulated polarization images. In our first 
paper only ADB was addressed; edge detection on the individual polarization 
images as well as DAB strategy were not considered. 

The results from our current study show that adapting edge detection meth¬ 
ods from Computational Intelligence to use in radar imagery is a new field 
worthy of exploration. In particular, our experiments show that modifying the 
gravitational method with Fu’s 9x9 neighbourhood produces better results 
than the unmodified method. They also show the importance of filtering when 
adapting edge detection techniques from optical to radar images. 

2. Basic concepts on SAR images 

Optical and SAR sensors measure the amount of energy reflected by a tar¬ 
get in various bands of the electromagnetic spectrum. The bands employed 
in most imaging radars use frequencies in the 2 MHz to 12.5 GHz range, with 
wavelengths ranging from 2.4 cm to 1 m. In this study, we used only the L-band 
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with wavelengths of [30 cm, 1 m] and frequencies of [1MHz, 2 GHz]. 

SAR systems generate the image of a target area by moving along a usually 
linear trajectory, and transmitting pulses in lateral looks towards the ground, 
in either horizontal (H) or vertical (V) polarizations [16], respectively denoted 
as H and V (see Figure [l]). In the past, the reception of the transmitted energy 
was made solely on the same polarization of the transmission, generating images 
in the HH and VV polarizations. Currently, with the advent of polarized and 
fully polarimetric radars (PolSAR - Polarimetric Synthetic Aperture Radar ), 
information about intensity and phase of the cross signals are also obtained, 
generating images relating to HV and VH polarizations. Usually, applications 
only consider the HH, VV, and HV polarizations. 




Figure 1: Horizontal and vertical signal polarizations transmitted by an antenna. Source: na 

The imaging can be obtained by gathering all the intensity and phase in¬ 
formation data from the electromagnetic signal after it has been backscattered 
by the target in a given polarization [18 . Each polarization in a given a scene 
generates a complex image, which can be thought of as two images, containing 
the real and imaginary values for the pixels, respectively. 

We denote the complex images from HH, VV, and HV polarizations as Shh, 
S H v, and Syv- Multiplying the vector [Shh Shv Syv] by its transposed 
conjugated vector [S^ H S^ v Sy V ] t , we obtain a 3 x 3 covariance matrix. The 
main diagonal contains intensity values; taking their square root, we obtain 
amplitude values. We denote the intensity images by Irh , Ihv, and Iyy and 
their corresponding amplitude counterparts by Ahh , Ahv , and Ayy. In this 
paper, we only considered the amplitude images, such as those depicted in 
Figure [2] 
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Ahh Ahv Aw 

Figure 2: Amplitude images for polarizations HH, VV, and HV from the same scene 


Speckle noise is multiplicative, non-Gaussian, and is proportional to the 
intensity of the received signal. Speckle degrades the visual quality of the dis¬ 
played image by sudden variations in image intensity with a salt and pepper 
pattern, as can be seen in Figure [2] It can be reduced with multiple looks in 
the generation of the complex images, causing degradation in spatial resolution. 
Another way to reduce noise is to employ filters, as will be discussed in the next 
section. 

In SAR image classification, one often uses samples from the classes in order 
to estimate the parameters of the distribution believed to underlie each class. 
Synthetic images can then be created using Monte Carlo simulation by taking 
the realization of the random variable associated to the class of each classified 
pixel. This artifice is useful to choose the most apt classifier for a given appli¬ 
cation: instead of relying solely on the original image, one takes the classifier 
that obtains the best average accuracy on the set of synthetic images. This 
methodology can also be used in other tasks, such as edge detection. 

3. Related Work 

One of the most successful edge detection algorithms for optical images was 
proposed by Canny [2], based on the following guidelines: i) the algorithm 
should mark as many real edges in the image as possible; ii) the marked edges 
should be as close as possible to the edge in the real image; iii) a given edge in 
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the image should only be marked once; and iv) image noise should not create 
false edges. It makes use of numerical optimization to derive optimal operators 
for ridge and roof edges. The usual implementation of this method uses a 3 x 3 
neighbourhood. 

A more recent multi-scale edge detection method was proposed by Lopez- 
Molina et al [8], using Sobel operators for edge extraction and the concept 
of Gaussian scale-space. More specifically, the Sobel edge detection method 
is applied on increasingly smoother versions of the image. Then, the edges 
which appear on different scales are combined by performing coarse-to-fine edge 
tracking. 

The gravitational edge detection approach was first proposed by Sun et al |3| 
and applied to optical images. It is based on Newton’s Universal Law of Gravity, 
described by Equation 0: 


h, 2 =Gx 


m i x m 2 



(i) 


where mi and m 2 are the masses of two bodies; r is the vector connecting them; 
/i } 2 is the gravitational force between them; ||.|| denotes the magnitude of a 
vector; and G is the gravitational constant. In the analogy proposed by Sun et 
al |3|; the bodies are the gray level values of pixels in a grid; G is a function 
of the values of the pixels in a given window; the distance between any two 
adjacent pixels is equal to 1; and, when computing the resulting force of the 
pixel in the center of a window; the pixels outside that window are considered 
negligible. Lopez-Molina et al [4 extended this technique, proposing the use 
of a Triangular Norm [19] in place of the product between the two masse^j by 
first normalizing the gray level values to [0,1]. The authors treat edges as fuzzy 
sets for which membership degrees are extracted from the resulting gravitational 
force on each pixel. They take G as a normalization constant, calculated so as 
to guarantee that the resulting forces lie in [0,1]. Also, in the normalization of 


1 Triangular norm operators are mappings from [0, l] 1 2 to [0,1], that are commutative, as¬ 

sociative, monotonic, and have 1 as neutral element. 
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gray level values into [0,1], a small value Sq is added beforehand to both the 
numerator and denominator so as avoid pixels with value 0, which would have 
too strong an effect on neighbouring pixels. The authors used 3x3 and 5x5 
windows as well as several prototypical triangular norms. 

The so-called Lee (or sigma) filter introduced in 1983 [20], is still in use today 
due to its simplicity, its effectiveness in speckle reduction, and its computational 
efficiency. It is based on the fact that, under Gaussian distribution, approxi¬ 
mately 95.5% of the probability is concentrated within two standard deviations 
from the mean. The filter estimates the mean and the standard deviation of 
samples around each pixel, and only those values within this interval are used 
to compute the local mean. Lopes et al ED proposed an adaptive version for 
this filter, here referred to as “Enhanced Lee”. 

Torres et al [12] recently proposed a nonlocal means approach for PolSAR im¬ 
age speckle reduction based on stochastic distances; the method can be tailored 
to any distribution, both univariate of multi-variate. It consists of comparing 
the distributions which describe the central observation for each pixel, and each 
of the observations which comprise a search region. The comparison is made 
through a goodness-of-fit test, and the p -value of the test statistic is used to 
define the convolution matrix which will define the filter: the higher the p -value 
the larger the confidence and, thus, the importance, each observation will have 
in the convolution. In Torres et al’s proposal, the tests are derived from h-<fi 
divergences between multi-look scaled complex Wishart distributions for fully 
PolSAR data [21]. Their results are competitive with classical and advanced 
polarimetric filters, with respect to usual quantitative measures of quality. 

Fu et al m proposed a statistical edge detector suitable for SAR images 
which uses the squared successive difference of averages to estimate the edge 
strength from the sliding window. An interesting feature of this paper is the 
proposal of a specific type of 9 x 9 neighbourhood, shown in Figure [3] In a pre¬ 
vious paper EE we proposed a modification of the gravitational approach using 
Fu et al’s neighbourhood: given a central pixel in a 3 x 3 window in an image, 
the values considered for the surrounding pixels in the window are no longer the 


ones in the original image, but the mean values in this new configuration. 



Figure 3: Standard 3x3 and Fu’s neighbourhood w 


4. Materials and Methods 

We compare the edge detection methods proposed by Canny [2] and by 
Lopez-Molina et al [8] to the modified gravitational approach using the prod¬ 
uct T-norm, followed by thresholding. The effect of preprocessing the images 
through filtering is also studied, using the filter described by Torres et al [12] and 
the Enhanced Lee filter m- We study the behaviour of Lopez-Molina’s method 
with the usual 3x3 window as well as a modified version of this approach, pro¬ 
posed in a previous paper m , involving the neighbourhood proposed by Fu et 
al [10]. 

The input for Canny’s and Lopez-Molina’s edge detector are images in, re¬ 
spectively, {0,..., 255} and [0,1]. Image values are, thus, mapped into these 
sets prior to edge detection. For the Lopez-Molina methods (the original and 
modified versions), we normalize further to [0,1], using Sq = 1 and making 
q' = (q + 1)(255 + l) -1 , where q and q' are the old and new value of a given 
pixel, respectively. 

4.1. Working image 

We apply the methods on data derived from a fully polarimetric image, pre¬ 
sented by Barreto et al [13], from an agricultural area in the Amazon region 
in Brazil (see Figure [T]). The authors describe a classification experiment using 
classes of interest from that area, such as water and different types of crops and 
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natural vegetation, at different stages of growth. Samples from the classes from 
band L are used to estimate the parameters of the complex Wishart distribution 
associated to each class. The results are assessed using a mosaic with the classes 
that was created using the derived Wishart distributions. Figure [5] illustrates 
the approach. For our study we apply the edge detection methods on twenty 
independently simulated mosaics amplitude images, using the parameters esti¬ 
mated in [T3] to assess the quality of the methods. 



Figure 4: Images derived from a scene in Bebedouro in Brazil (not registered): a) Landsat 
RGB composition and b) SAR L-band RGB composition (source: |13|4 



Figure 5: Images derived from a scene in Bebedouro in Brazil: a) training samples used to 
generate Wishart distributions and b) synthetic mosaic images generated using the Wishart 
distributions estimated in m from image samples (source: m) 
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4-2. Quality assessment 

The quality of the results is assessed by the Baddeley’s Delta Metric (BDM) 
0 , by comparison with what would be the perfect result, discarding those 
pixels close to the outer frame. 

Let x and y be two binary images, seen as mappings from A to [0,1], where 
A is a set of sites arranged in a grid (positions). Let p be a metric on A, such 
as the Euclidean distance, and d(i, A) be the distance between a site i and a set 
A C A, defined as 

d(i, A) = minp(i,jf). 
jeA 

Let 6(x) = {i G A | X} = 1} denote the set of foreground sites in x. BDM 
between x and y, denoted as A PjWm (., .), is then defined as 

A p , w (x, y) = L ^ | w(d(i, fo(x)) - w(d(i, b( y)) \ p j , 1 < p < oo (2) 

where w is a strictly increasing concave function satisfying u>(0) = 0. Here we 
use w(t) = t and p = 2, as in [4]. 

Throughout the text, we display BDM results in [0,100] instead of [0,1], for 
the sake of readability. 

5. Proposed methodologies 

Edge detectors use a window around a center pixel to verify whether that 
pixel belongs to an edge or not. When adapting optical image edge detectors 
to radar imagery, we have to find the means to deal with speckle. The main 
contribution of this study is to modify the original 3x3 window used by the edge 
detection method proposed in [4 for use in radar imagery such that the value in 
each cell in the window is no longer the original one but the aggregation of set 
of neighbouring pixels, according to a larger 9x9 non-standard neighbourhood 
proposed by Fu et al m- We here investigate this particular combination of 
method and neighbourhood, but the same procedure can be applied using other 
edge detection methods and/or non-standard filters. 
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Frequently, a single band is used in edge-detection, resulting in a gray level- 
image that is then binarized at some point (the usual implementation of some 
methods, like Canny’s, already involve a binarization step). In radar imagery, 
very often one deals with more than one band at the same time (e.g. intensity 
images coming from different polarizations, or complex images in the fully po- 
larimetric case), aiming at using the richness of information to compensate for 
the speckle noise. Therefore, the question of when to aggregate results has to 
be addressed. One may, for instance, first aggregate the bands and then apply 
the edge detector on the aggregated image, or else apply the edge detector on 
the individual bands and then aggregate the edge images. These two methods 
usually yield different results. 

Here, we propose a typology for experiments using radar imagery, consid¬ 
ering different orderings of three steps: edge detection on gray level images, 
binarization of gray level images, and aggregation of results. In the aggregation 
step, the input may be either gray level or binary images, depending on whether 
the binarization is made immediately after edge detection or not. Three strate¬ 
gies can then be envisaged to perform edge detection experiments with radar 
images: 

• DAB (edge Detection on non-binary images, Aggregation of the resulting 
non-binary images, Binarization) and 

• ADB (Aggregation of non-binary images, edge Detection on the resulting 
non-binary image, Binarization). 

• DBA (edge Detection on non-binary images, Binarization, Aggregation of 
the resulting binary images). 

Options ABD, BAD and BDA are not considered, since that would mean ap¬ 
plying edge detectors on the binary images. 

In this work, we focus on the DAB and ADB strategies. For both of them, 
we use the arithmetic mean to aggregate gray level images. Strategy DBA, 
involving the aggregation of binary images, is left for future study. 
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When no aggregation is considered, the strategies are reduced to only edge 
detection and binarization. For example, for the HH, HV and VV polarizations, 
we obtain strategies DB-HH, DB-HV and DB-HH. Note that some methods 
already incorporate the binarization step in the edge detector. That is for 
instance the case of all methods discussed previously. However, to be consistent 
with the notation, we will denote by ADB the strategy in a method that includes 
binarization, such as Canny’s and the multi-scale method, when it is applied to 
the gray level image resulting from the aggregation of the images from the HH, 
HV and VV polarizations. 

6. Experimental Results 

The output of the Lopez-Molina gravitational method is an image with val¬ 
ues in [0,1]. In order to obtain binary indicators of edges, the authors use a 
hysteresis transformation. Here, we use a simple threshold and search for the 
value in the [0.05,0.15] interval which produces the best BDM. For Canny’s 
method, we search for the best value for the noise standard deviation parameter 
a in the interval [0.3,1.5]. The intervals above for both Canny and Lopez-Molina 
are the ones that presented the best results by trial-and-error. The following 
parameters were used for the Lopez-Molina multi-scale method, as suggested 
in 0: 5 a = 0.25, cr e {0.50, 0.75,1.00, 1.25,1.50,1.75, 2.00, 2.25, 2.50, 2.75, 3.00, 
3.25,3.50,3.75,4.00}. 

We applied two filters (Torres [12| and Enhanced Lee [Tl]) on intensity val¬ 
ues, which were then transformed in amplitude before further processing. 


Table 1: Average BDM results for Canny’s method, with standard deviation inside parentheses 


Strategy 

No filter 

Torres filter 

Enh. Lee filter 

DB-HH 

26.72 (1.22) 

23.40 (1.92) 

28.85 (1.86) 

DB-HV 

30.70 (1.39) 

26.74 (1.33) 

66.40 (0.42) 

DB-VV 

29.81 (2.17) 

26.28 (1.09) 

24.83 (1.98) 

ADB 

27.87 (1.16) 

48.16 (0.003) 

30.24 (1.20) 
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Table 2: Average BDM results for the multi-scale method, with standard deviation inside 
parentheses 


Strategy 

No filter 

Torres filter 

Enh. Lee filter 

DB-HH 

28.23 (0.89) 

28.00 (0.93) 

25.36 (2.11) 

DB-HV 

28.34 (1.13) 

24.55 (2.97) 

19.56 (2.36) 

DB-VV 

25.62 (1.10) 

28.42 (0.56) 

28.89 (2.25) 

ADB 

25.15 (3.23) 

24.37 (1.52) 

20.71 (3.97) 


Table 3: Average BDM results for the gravitational method; with standard deviation inside 
parentheses 


Strategy 

No filter 

Torres filter 

Enh. Lee filter 

DB-HH 

33.89 

(1.98) 

26.61 

(2.00) 

38.97 

(0.79) 

DB-HV 

31.95 

(0.46) 

27.14 

(1.22) 

32.26 

(2.78) 

DB-VV 

32.35 

(1.47) 

28.95 

(0.95) 

43.65 

(1.13) 

DAB 

29.26 

(1.62) 

25.91 

(1.55) 

27.71 

(2.62) 

ADB 

31.50 

(0.82) 

26.63 

(1.27) 

18.24 

(3.41) 


Tables Si! and [4] show the results for BDM mean and standard devia¬ 
tion after applying four methods to twenty simulated mosaic images: Canny’s 
method, Lopez-Molina et al’s multi-scale method, Lopez-Molina et al’s original 
gravitational method, and Lopez-Molina et al’s method modified using Fu’s 9x9 
neighbourhood. 

We see that the best BDM average values were obtained with the use of 
Lopez-Molina et al’s gravitational method modified by Fu’s neighbourhood, us¬ 
ing the ADB and DAB strategies, both preprocessed with the Enhanced Lee 
filter. Both are significantly higher than the other procedures. 

In Table[l] we see that filtering did not have a significative impact on Canny’s 
detector. The same is true, to a lesser degree, for most results of the multi-scale 
and (unmodified) gravitational methods, as can be seen in Tables [ 2 ] and i In 
these methods, there is a slight advantage in preprocessing the images using the 
Enhanced Lee filter. However, filtering has an impressive effect on the Lopez- 
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Table 4: Average BDM results for the gravitational method modified by Fu’s neighbourhood, 
with standard deviation inside parentheses 


Strategy 

No filter 

Torres filter 

Enh. Lee filter 

DB-HH 

25.27 

(0.76) 

22.18 

(0.48) 

17.79 (3.05) 

DB-HV 

26.48 

( 1 . 00 ) 

24.21 

(0.65) 

18.40 (5.75) 

DB-VV 

21.41 

(1.97) 

18.14 

(0.77) 

17.83 (2.54) 

DAB 

22.67 

(2.29) 

18.97 

(1.62) 

5.43 (1.68) 

ADB 

23.80 

(2.23) 

22.74 

(0.50) 

5.16 (0.36) 


Molina gravitational method modified with Fu’s neighbourhood. In particular, 
the best results are obtained for strategies DAB and ADB with preprocessing 
with the Enhanced Lee filter. 

Figure [ 6 ] shows the negative images corresponding to the best results, accord¬ 
ing to BDM, obtained by the edge detection methods and the filtering strategies 
with the best average values; note that the image boundaries are depicted only 
for illustrative purposes. We see that according to BDM the best binary image 
(depicted in Figure [ 6 Jd) presents little noise and most of the regions are sepa¬ 
rated, even though the lines are rather thick. We also see that BDM was able 
to distinguish the best image from the others. 

Figure [ 6)3 shows the best results from the methods come from filtered images, 
which raises the question of how important preprocessing by filtering is. In what 
follows we discuss the details of the gravitational method using the original 3 x 3 
and Fu’s 9x9 neighbourhood in relation to filtering. We take the simulation that 
obtained the best BDM results for each type of neighbourhood. We see in these 
examples, that filtering does, indeed, ameliorate the results for all methods. 

Figures[7]a),[7|b) and[7]c) respectively show that: using the 3x3 neighbour¬ 
hood for the original gravitational method, the unfiltered image is very noisy; 
the Torres filter reduced the noise and separated the regions; and the Lee filter 
detected false edges. In Figure [ 8 j we see that Fu’s 9x9 neighbourhood detected 
almost all the edges, especially when using Lee’s filter. Filtering for the modified 
method presented a larger trade-off between detection of edges and reduction 
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Figure 6: Best BDM results obtained from the best methods (average): a) Canny (DB-HH, 
with Torres filtering, BDM = 18.51), b) Multi-scale (DB-HV, with Enh. Lee filtering, BDM 
= 14.94), c) Gravitational (ADB with Enh. Lee filtering, BDM = 10.96) and d) Gravitational 
and Fu (ADB with Enh. Lee filtering, BDM = 3.05) 


of noise (some edges were detected using Torres filter with an increase of noise 
when compared to the unfiltered image). 

When we compare the results in Figures [7] and [8] we see that the gravitational 
method modified with Fu’s 3x3 neighbourhood clearly produced better results 
than the method with the original 3x3 window, which agrees with the BDM 
evaluation. 

Figures [9j [lOj and [TT] show the best edge detection results obtained for 
the Bebedouro SAR image in terms of visual analysis, with the application 
of the Enhanced Lee filter, for the parameterizations used here. The figures 
present the results for HV polarization; in general, HH (respec. VV) polarization 
produced images with more noise (respec. less information) than HV. Figures [To| 
and EH respectively, show the results of the application of ADB and DAB 
aggregation strategies. ADB in general produced binary images with very little 
information for all methods; the decrease of noise in relation to the individual 
polarizations does not compensate the lack of information. In general, the 
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Figure 7: Results for the gravitational method with original 3x3 window, on a single simulated 
image from ADB: a) no filtering (BDM = 31.90), b) Torres (BDM = 27.73) and c) Enh. Lee 
(BDM = 10.96) 



Figure 8: Results for the gravitational method modified with Fu’s neighbourhood, on a single 
simulated image from ADB: a) no filtering (BDM = 23.80), b) Torres (BDM = 23.76) and c) 
Enh. Lee (BDM = 3.05) 

DBA aggregation method produced results with less noise for the gravitational 
method, with and without modification, than the results obtained with the 
individual polarizations. The best results for the gravitational method, both 
with and without modification with Fu’s neighbourhood, were obtained with 
thresholds around the same interval that produced the best results using the 
mosaics. 

7. Conclusions 

Contrary to what happens with optical imagery, few algorithms are specif¬ 
ically dedicated to PolSAR image edge detection [10,. One interesting means 
to create edge detection algorithms for SAR images is to modify those created 
for optical images, in such a way as to reduce the non-Gaussian noise. Here 
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Figure 9: HV Bebedouro binary images, with Enh. Lee filter: the first row depicts results 
obtained using the Canny and the Multi-scale methods; and the second row depicts results ob¬ 
tained using the original Gravitational method and its modification with Fu’s neighbourhood, 
(the latter two methods use binarization threshold=.2) 

we have investigated the modification of a method issued from Computational 
Intelligence for optical imagery, the gravitational edge detection method exten¬ 
sion proposed in [4] (see also 0 ), to synthetic aperture radar imagery. In order 
to deal with speckle, we modified the gravitational method with a non-standard 
9x9 neighbourhood configuration proposed by Fu et al m- considering a 
3x3 window centered around a given pixel, the value of any pixel in the win¬ 
dow becomes the average value of the region associated to that pixel in the 
non-standard neighbourhood configuration. 

Considering that SAR imagery has different polarizations, and that their 
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Figure 10: ADB Bebedouro binary images, with Enh. Lee filter: the first row depicts results 
obtained using the Canny and the Multi-scale methods; and the second row depicts results ob¬ 
tained using the original Gravitational method and its modification with Fu’s neighbourhood, 
(the latter two methods use binarization threshold=.2) 

joint use may compensate for the presence of speckle, we also proposed a ty¬ 
pology of experiments regarding aggregation of these images. In particular, we 
addressed two procedures: DAB (edge Detection on non-binary images, Aggre¬ 
gation of the resulting non-binary images, Binarization) and ADB (Aggregation 
of non-binary images, edge Detection on the resulting non-binary image, Bina¬ 
rization). 

For means of comparison, we also addressed the use of two other edge de¬ 
tector methods stemming from the realm of optical images: the traditional 
method proposed by Canny [2] and a recent multi-scale one coming from Com- 
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Figure 11: DAB Bebedouro binary images, with Enh. Lee filter: the first and second rows 
respectively depict results using the original Gravitational and Gravitational modified with 
Fu’s neighbourhood, the first and second columns respectively depict results obtained with 
binarization thresholds .1 and .2 


putational Intelligence, based on Sobel operators for edge extraction and the 
concept of Gaussian scale-space [8]. 

We studied the effect of filtering the images prior to edge detection by two 
procedures: Enhanced Lee m and Torres [12] filters. The methods were applied 
on twenty samples of a scene, which were simulated using Wishart distributions 
derived from a fully polarimetric image m- Using both visual inspection and 
the Baddeley Delta metric im we verified that the combination with the Lopez- 
Molina technique with the 9x9 neighbourhood proposed by Fu et al m and 
preprocessing with the Enhanced Lee filter produced the best results. 
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This paper is an extended version of m ; together, these studies represent 
a first step towards investigating the use of edge detection methods derived 
from Computational Intelligence techniques for use in SAR images. The main 
implication of our results is that the joint use filtering and neighbourhood mod¬ 
ification on those methods, as well as the use of aggregation of the individual 
polarization images, are able to deal with speckle, which is crucial when detect¬ 
ing edges in radar imagery. 

Future work includes modifying the Lopez-Molina method with other types 
of neighbourhoods, such as Nagao-Matsuyama [22]; to verify the performance 
of other T-norms than the product to calculate the gravitational forces; and to 
perform preprocessing with other filters. We also intend to investigate the use 
of the proposed procedure with other edge detection methods, such as the one 
described in [23] . involving fuzzy sets. 

We would like to better address the issue of aggregation. Here we have dealt 
exclusively with the aggregation of non-binary images, using the arithmetic 
means in strategies DAB and ADB. In the future, we intend to explore aggrega¬ 
tion of the images considering families of operators in general, such as weighted 
means, ordered weighted means (OWA), T-norms, and T-conorms [T9] . Also, 
we intend to study other operators than the average to perform aggregation of 
pixel values in regions of a non-standard neighbourhood. 

Moreover, we intend to assess the results using other methods than BDM, 
such as the one proposed recently by Frery et al m- We would also like to draw 
comparisons with other edge detection algorithms, such as the one proposed by 
Fu in 2012 [10]. 

Last but not least, we intend to verify the use of the approach consider¬ 
ing fully polarimetric images (PolSAR), instead of just intensity images. In 
this case, Torres filter, designed specifically for PolSAR images, can be more 
adequately employed. 
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